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Abstract 

We study the problem of spin diffusion in magnetic systems without long- 
range order. We discuss the example of the ID spin chain. For the system 
described by the Heisenberg Hamiltonian we show that there are no diffusive 
excitations. However, the addition of an arbitrarily small dissipation term, 
such as the spin-phonon interaction leads to diffusive excitations in the long 
time limit. For those excitations we estimate the spin-diffusion coefficient by 
means of the renormalisation group analysis. 
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I. INTRODUCTION 



The spin dynamics in magnetic systems without long-range order is a long-standing prob- 
lem. It is generally believed that in the high temperature limit, where no long-range 
order is present, the microscopic spin fluctuations are governed by the diffusion equation, 
i.e. that at small frequencies and momenta the spin correlation function has a diffusive pole 



where D is the diffusion constant and Z is the residue. 

Although this idea of spin diffusion is quite common, we are not aware of any theoretical 
approach, within which one could "derive" the correlation function Eq. (|I]), starting from 
the nearest neighbour Heisenberg Hamiltonian 



In this paper we are presenting such an approach for ID spin-1/2 Heisencberg chains 
of the infinite length. At any non-zero temperature the chain is in disordered state with 
exponentially decaying spin-spin correlations and even at T = the correlations decay as 
a power law, so there is no true long-range order. We show that in one dimension it is 
impossible to derive the diffusive form of the spin-spin correlation function Eq. (|l]) from the 
Heisenberg Hamiltonian Eq. (|^), without any kind of additional dissipation mechanism. We 
also present a general argument supporting this result. We show then, that if one takes 
into account an additional dissipation, for example due to spin-phonon interaction (which 
is always present in any real system at finite temperature), then the renormalisation group 
approach leads to the correlation function Eq. (|l|) in the long-time asymptotic. 




(1) 
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Our results could be applied to the materials like KCuFs, CUSO4 ■ 5H2O, Sr2CuOz. In a 
broad temperature range they are nearly ideal ID antiferromagnetic chains with the coupling 
constant J ranging from 1.45 K in CuSO^ ■ 5H2O |1[ to 190 K in KCuFs 0§ and 1300 
K in the novel Sr2CuOz 0- The diffusion equation was successfully used in a number of 
experimental papers to fit the data and explain the results of the experiments This 
approach was also confirmed by computer simulations |^. 

The rest of this paper is organized as follows. In Section II we briefly review the map- 
ping of the spin system onto ID fermions. In the Section III we map the fermions onto 
a boson system, justifying bosonisation procedure at finite temperatures. Section IV gives 
the calculation, which shows the absence of spin diffusion in the Heisenberg model Eq. 
and presents some general qualitative arguments, supporting this conclusion, based on the 
theory of the Sine-Gordon equation. Section V discusses the spin-phonon interaction and 
its renormalisation. Section VI is conclusions and discussion of results. The details of the 
calculations are presented in Appendices. 



II. FERMION MODEL 



In this section we review the well-known procedure of mapping the Heisenberg model 
Eq. (0) on a system of spinless fermions and establish the notations. 

The spin model Eq. (0) can be transformed into a model of spinless fermions, noting that 
operators Sf and anticommute. The Jordan- Wigner transformation then relates spin to 
fermion operators {a^,a\) via 

Sf = a\exp {^-K ^ a]a^ , 

St = aU, - ^ (3) 
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When transforming the Hamiltonian Eq. (|]) the spin-flip terms give rise to the motion of 
the fermions (kinetic energy in the fermion Hamiltonian) and ■ S^j^^ interaction leads 
to a fermion-fermion interaction between adjacent cites. Since 5*^ is quadratic in fermion 
operators, the interaction between fermions is the four-particle interaction. Since the original 
spin model was formulated on a lattice, all possible types of four-fermion interaction are 
present, including the umklapp term. 

The ID fermion models are often treated using bosonisation In the case of massless 
fermions with the four-fermion interaction with small momentum transfer it allows exact 
solution. It is shown, that in the thermodynamic limit that system can be mapped onto a 
system of free bosons. The propagator of a free boson has a pole aX uj = ck and therefore 
corresponds to the particle, propagating without dissipation, so that this interaction can 
not lead to any signiflcant change in the long-time dynamics of the system. Therefore the 
only source of dissipation might be interaction with large momentum transfer, such as the 
umklapp term. 

The treatment of the umklapp term is much more complicated because in the boson 
language it corresponds to the highly non-linear term, namely cos(2/30) (where 0(x, t) is the 
boson fleld and the number /? is a parameter of the transformation, depending on the fermion 
interaction with small momentum transfer). This term does not conserve momentum and 
is only marginally irrelevant, so at T > it might lead to some new dynamics. Thus 
we should investigate the impact of that term on the long-time fermion (and thus spin) 
dynamics, disregarding all other possible four-fermion interaction terms. The next two 
Sections present the results of this investigation. 

Since we are interested only in the low-energy behaviour of the system, we can linearize 
the fermion kinetic energy, thus allowing exact bosonisation. Therefore our model Hamilto- 
nian is 



(4) 



where is the operator of a "right" ("left") mover. The quantity we are interested in 

here is the density- density correlation function, which corresponds to the {S^ ■ Sj) correlator 
in the spin problem. We shall now investigate, whether the long-time asymptotic of this 
correlator has the diffusive form, Eq. (|1|). 



III. BOSONISATION 



As we mentioned above our main technical tool will be bosonisation. The procedure 
is well established in ID at zero temperature. We consider finite temperatures and we are 
looking for the long-time asymptotic of the spin-spin correlator. It is quite difficult to obtain 
the results in that limit in the Matsubara technique, since the analytic continuation from 
the Matsubara frequencies u;„ to real frequencies much less that the inverse temperature 
would require a precise knowledge of the Green's functions on the infinite range of Un, which 
is usually not the case in perturbation theory. Therefore we have to resort to the Keldysh 
technique [jlO|, which incorporates finite temperatures and real-time representation. In this 
Section we construct the bosonisation procedure for the Keldysh technique and then confirm 
it's correctness by comparing the results of the perturbation theory for bosons and fermions. 
The fermion and boson Green's functions in space-time representation, which we are using 
in our calculations, are presented in Appendix A. 

Thirty years ago Keldysh has presented a field theoretical technique to calculate the 
real-time correlation functions of a quantum system. To allow the treatment of advanced 
and retarded correlators, Keldysh introduced the time contour C (Fig. 1) with the upper 
branch going in positive direction from — oo to +00 and the backward lower branch. All 
the operator products are now time-ordered along the contour C. To distinguish particles 
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on the upper and the lower branches of the contour the fermion field operator if)^ is given 
an index 7, which equals to 1 on the upper branch and 2 on the lower. Green's functions 
become 2x2 matrices with respect that index. 

In a one dimensional fermion problem we have four different operators ■j/'^^^^ - left and 
right movers on both branches of the contour. Since operators on each separate branch 
are completely analogues to the zero-temperature operators, whe can proceed with the 
bosonisation separately on each branch in exactly the same way as at T = 0. Thus we 
introduce two boson fields, 0^, (one for each branch), which we shall treat as two components 
of the Keldysh field. The resulting bosonised Hamiltonian will thus be formulated in the 
Keldysh technique also. 

The left 0^ and right (p^ moving bose fields expressed via (p^ and its canonically conjugate 



(x) =F I P^{x')dx' 
J —00 



•^0 2^T^J2\^\ 



,(p)e*^'" + h.c. 



(5) 



As in the usual procedure 0;^^^^ are functions of only (x^t). 

The fermion operators are constructed in analogy with the zero temperature case 




where = Att and the upper sign corresponds to the left mover. 

The commutation relations between fermion fields hold by exactly the same reason ||^ as 
at T = 0. The fact that we have a different time contour (the Keldysh contour C as opposed 
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to the usual time axis) does not change the calculation, for the integrals involved in Eq. 
are over space coordinate and the bose fields commute no matter which branch of the time 
contour they are on. Another way of saying this is that the Keldysh operators on different 
branches still correspond to the same particle. Dividing the time contour in two parts is a 
matter of mathematical convenience rather than physical distinction. 

The cut-off a is a lattice spacing and so should be the same for both bosons and fermions. 
The operator equality Eq. means that any correlation function (in the limit a — 0), 
calculated in the fermi theory with the cut-off a, is reproduced in the bosonic theory with the 
same cut-off if the fermion operator in the left-hand side is replaced by the bosonic operator 
in the right-hand side of Eq. (|^) . Using this operator equivalence we can construct the boson 
Hamiltonian from the fermion Hamiltonian Eq. (Q) . The fermion kinetic energy corresponds 
to that of bosons. The umklapp interaction term gives rise to the cosine interaction of the 
boson field. The conjugate operators P+ and P_ cancel out, so the interaction is a function 
only of the boson fields 0+ and 0_ itself, as it is in the zero-temperature bosonisation. The 
interaction constant V now acquires the factor from the prefactor in Eq. (^. The boson 
Hamiltonian therefore is 



HB = {d^(pf + V'cos2l3(j). (7) 



where V = 



a 



Due to the special form of the interaction (cos 2/30) it is most convenient to formulate 



the Keldysh technique in the path integral representation, developed by A. Schmid |jrT[ . The 
boson action for our model is 



A = dxdt 



0!R>,- \/'(cos2/?0_,-cos2/50_) , (8) 
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where 0^ and 0_ are the two components of the Keldysh boson field, 7 and v denote Keldysh 
indices. B^^ is the kinetic operator, so that the inverse is proportional to the boson 
Green's function in the Keldysh basis 



K' = -^\ ^ 1 (9) 



We shall later use these functions in the space-time representation (see Appendix A). 

The fermion density in the conventional bosonisation is represented by the spatial deriva- 
tive of the boson field 



P = (10) 



Here in the same manner we write the fermion density operators. The density-density 
correlators are represented by functional integrals in which the predexponential is some linear 
combination of the spatial derivatives of the bose fields corresponding to that particular 
correlator. The advanced correlator is 

{pp)a = n{xi - X2) = ^ J V [4>i,^2] di(f)l{xi)T^pdi4)f^{x2) exp(-^), (11) 
where the density vertex TaB is 



1 



1 -1 



Ta^ = - . (12) 



TT 



1 -1 



We now check, following the approach of R. Shankar [Q, that the above construction 
yields correct results in perturbation theory. The density- density correlation function for 

8 



non- interacting fermions should be calculated separately for the "left" and "right" movers 
and the results should be added. That gives for the advanced correlator 



2 

no = - , \ , (13) 



where 8 is infinitesimately small. In the next order we have two topologically nontrivial 
diagrams given on Fig. 2 Again we have to repeat the calculation for the "left" movers 
and add the results with proper combinatorial factors. For simplicity we give here only the 
imaginary part of the first-order correction, which for us is most important 



Mn, . |1 ^ (14) 



The next order of perturbation is discussed in the next Section. 

We now turn to the calculation of the same perturbation series on the boson language. In 
the non-interacting case we have to calculate a simple Gaussian integral which in momentum 
space gives 



i7o(u;, q) = -q' Da{uj, q) = - — , (15) 

vr T[ uo^ — q'^ — liouo 



which is the same as Eq. (|TB|). DA^oj.q) is the advanced boson Green's function, connected 
to the original Keldysh basis via 
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Da = D— - D+- 



(16) 



We now expand the exponent in the integral Eq. ([TT|) in the series in V. In the first 
non-trivial order we get 



exp J dxdt 0* 



cos2/30+(yi)-cos2/50_(yi) cos 2/30+ (7/2 )- cos 2/30_ (2/2) • (17) 



Since the cosines contribute the linear terms in the exponent, the integral is still Gaussian 
with the same prefactor of the exponent. This feature will remain in all higher orders of the 
perturbation series. The calculation of that integral we shall discuss in more detail due to its 
importance for the later arguments. In the first-order integral Eq. (0) we have four different 
terms of the same type, arising from the cosines. In the exponent they have combinations 
like 2/3(0^(^2) ^02(2/1)) with all possible permutations of indices. To calculate the functional 
integral we make the Fourier transform of the bose fields. We perform the transform in the 
most general way, since the same expressions will appear later. As we show below in any 
order of the perturbation theory we shall need to calculate averages of the form 



(9i0;(a;i)r„^ai0^(x2) ^xp (^/ ^ ^""^r^iKVi) <P,{k) + h.c. j). (18) 



For the first-order integral Eq. ([T7|) we four different factors ^^^I*{k): 
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(19) 



Completing the square in the exponent we calculate the functional integral in Eq. ([T7|) and 
apart from the numerical factor we get in the momentum space 

n,{uj,q) ~ J d'k kq J n A^E B^^^k) (^^/.(fc, (20) 

{r} 



where 



K,{m) = (/ W/;(A;',y,)5-/(fc') W/|(A;',y,) j , (21) 



and the sum is over four different factors ^^^Ij^{k,yi) corresponding to the four different 
terms of the type (|TBp in the first order integral Eq. (PTj). 

Here the integration measure Yid'^Vi is equal dyidtidy2dt2. All dependence on these 
variables is contained in the factors I*{k,yi). These factors are nothing but exponentials 
^gjfcy -^vith different signs. Therefore the integral in the exponent in Eq. (|2T|) is just the 
Fourier transform of the boson Green's function back from momentum space to the real 
space. The boson functions in real space are logarithms, so all four exponents Ki{xi) can 
be evaluated. All four have the same structure 

where 
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S'\y^ t) = —T-2 — . ^""^""^ . 2 — ; (23) 

sinn TTT[y — t) sinh TcT[y + t) 



while f{y,t) is some algebraic function, which has no poles. In the long-time asymptotic 
Ki{y,t) acts like a derivative of a delta-function. It has (due to the denominator Eq. ([23|)) 
a sharp singularity at y,t = and decays very fast as the variables go away from that 
point. Therefore the integral over y,t is dominated by the small region of the size a around 
the origin. The remaining spatial integration over the sum yi + y2 involves only functions 
^^^I*{k, yi) from the prefactor since Ki{yi) depends only on the difference of the variables yi 
and 2/2- Since (/c, yi) are exponentials the integral is easily evaluated to yield 5{k — q), 
which solves the momentum integration. So after the summation over all four terms we get 



the result Eq. (0), as expected. 

The next orders of the perturbation series, which we discuss in the following Section, 
could be calculated in the same manner and are also the same for the boson and the fermion 
formulations of he theory. The bosonisation procedure is thus justified. 



IV. ABSENCE OF DIFFUSION FOR THE HEISENBERG MODEL 



Although the boson and the fermion versions of the theory are completely equivalent and 
should give the same results in the perturbation theory, the bosonised version allows easier 
evaluation or the higher orders of the perturbation theory. Indeed, all the higher order cor- 
rections differ from the first order Eq. (p^ only by appearance of additional cosine brackets 
(and corresponding spatial integrations). That means that in any order of perturbation we 
have to calculate averages of the form Eq. (|18|). The functions ^^^If^{k, Xi) will be now sums 
of exponentials, for example in the next order of perturbation we will have terms like 
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(24) 



similar to the first order factors Eq. ([T9|), but constructed from four different exponentials. 
We can still perform the functional integration and get the same general formula Eq. (pO|). 
Since functions ^^^I^(k,Xi) are exponentials the integral in the exponent in Eq. (|20|) still 
yields the boson Green's functions in the space-time representation, but now instead of one 
such a function as in the first order we have a sum of them. In the second order we shall get 



where S{x) is the "singular" denominator, defined in Eq. (p^). 

In higher orders functions K have the same structure, but with more factors S{x) in the 
numerator and denominator. These factors have the same singular behaviour as described 
in the previous Section. Therefore the integration over the multidimensional space will be 
dominated by the regions where two pairs of variables are almost equal to each other, for 
example Xi = X2 and X3 = X4. The integration over one such pair effectively reduces the 
function Eq. ( p5D to the form of the previous order Eq. (p^ (in the long-time asymptotic). 
One can see it by inspecting Eq. (pSj). Consider, for instance, the contribution from the 
region xi = X2. In this region K2 ~ S~^{xi — X2)5'~^(x3 — X4), so the integral over X2 leads 
to the same form (~ 1/ S{x^ — 0:4)), which is the main part of the first order function Ki 
Eq. (|22|). The integration over X3 and x^ is then the same as in the first order . Therefore 
the spatial integration in in the higher order terms does not yield any additional singularity 
to that produced by the prefactor of the exponent, which is the same as in the first order 
Eq. (^) and contains just two boson Green's functions. So the structure of any term in the 
perturbation series is 



K2{xi, ...Xi) 



S{X2 - X3)S{xi - X4) 



/2(Xi, ...X4) 



(25) 



S{xi — X2)S{xi — X3)S{X2 — X4)S{X3 — X4) 
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Imn,{uj,q) ~ P giuj,q). (26) 

(u;^ — g^) 

Here g{uj,q) is some function, which has no poles at small uj,q. We note that Imni{u,q) 
differs from the first order Eq. ([T7|) only by g{uj, q), and P denotes principal value. 

Combining our conclusions we can write down the structure of the density-density cor- 
relation function in our model. 

TT uj^ — q^ + 2iduj {uj^ — q^y 

The significance of this result for us is in the fact, that the higher order corrections do not 
acquire additional poles in the imaginary part, which could sum up in a diffusive pole in 
n{uj,q). 

The fermion perturbation theory gives the same results in the low orders of the perturba- 
tion. In the second non-trivial order we have six topologically different diagrams, presented 
on Fig. 3. Our purpose is to show that the correction to the density-density correlation 
function, which is the sum of these diagrams, is not qualitatively different from the first 
order result. The only diagrams that produce the extra pole are the diagrams a and h on 
Fig. 3. They cancel each other exactly, so that the second order correction has the pole of 
the same order as the first order one. This cancellation seems to be an accidental property 
of the problem in the fermion representation, but the bosonisation approach shows that it 
happens in all orders of perturbation theory. 

Our results clearly show the absence of spin diffusion in the Heisenberg model. Instead 
of the diffusive pole in the density-density correlation function we found some kind of a 
propagating behaviour. Note that the second order pole in Eq. ( P7D should be regarded as a 
principal value, so the imaginary part of general susceptibilities does not contain unphysical 
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contributions proportional to — q^). This result should have been expected. In the 

boson representation our problem is essentially the Sine-Gordon model. It is known in the 
theory of the Sine-Gordon equation that due to the infinite number of conserved charges 
the excitations of the model are propagating and not diffusive. So by showing the mapping 
of the Heisenberg model onto the Sine-Gordon model we have showed the absence of spin 
diffusion in the model Eq. (|]). 



We now try to make it a bit more realistic. Even in the absence of disorder, at finite tem- 
peratures there always are phonons in the system. The exchange integral depends in general 
on the instantaneous separation of magnetic ions. We consider the linear approximation 
for spin-phonon coupling. That is we expand the separation-dependent Heisenberg coupling 
constant to the first order in ionic displacement. u{Ri) 



where Ri is the equilibrium position of the ion and ri = Ri + u{Ri). That leads to the 
simplest form of the spin-phonon interaction Hamiltonian 



V. SPIN-PHONON INTERACTION 



The Heisenberg model Eq. (0) does not exactly describes the physics of a real material. 



J(r, - r^) = Jo(i?i - Rj) + X HRi) - u{Rj)) , 



(28) 




(29) 



<ij> 
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We shall now try, treating interaction Eq. (^) as a small perturbation of the original 
Hamiltonian Eq. (0), to examine it's impact on the long-time spin dynamics. To do that we 
first bosonise Eq. (^). We get 



H,_,, = J'E|0,|' (bl + b,). (30) 



The boson self-energy in the first non-trivial order in Tit-ph is presented by the diagram on 
Fig. 4. The solid line represents a boson and the dashed - a phonon. For the imaginary 
part we are interested in we get 



ImT,{u!, q) 



riojavp , if q > T 

(31) 

r]u}T , if q < T 



where restoring the units t] ~ [J' /cf' and c is the speed of sound. The momentum depen- 
dence of ImL{uj^q) arises due to the momentum in the numerator of the phonon Green's 
function 



2 2 

-Dphiu^, q) = ^, (32) 



where S is infinitesimally small. Note that although we consider ID spin chain, the phonons 
in a real material are three-dimensional, so when evaluating the self-energy Eq. (|31|) one 
must integrate out the two other components of the phonon momentum. 

We shall now investigate, how the umklapp term renormalizes this self-energy. We now 
divide the boson field into "slow" and "fast" parts, separated by some cut-off ko. We 
integrate out all "fast" degrees of freedom (those with momentum bigger the cut-off) and 
see how the imaginary part of the boson self-energy (namely the coefficient 77) changes with 
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the cut-off. The result is that when the cut-off is bigger than the temperature, T, t] rises as 
some power of the cut-off. But after the cut-off becomes smaller than the temperature, the 
imaginary part grows exponentially. That means, that on the large distances the motion 
becomes diffusive. We notice that it should happen in any experimental situation, because 
it is just the existence of the spin-phonon interaction (no matter how small) is responsible 
for the diffusion. 

After we introduce the "slow" and "fast" variables as 



(33) 



The integral Eq. (0) then is 



n{xi-X2) = - I V[<Pl,<p^,]V 4,4 5iC(a^i)^a/35i0K^2) exp(-^) 



(34) 



where the prefactor contains only slow degrees of freedom since we are looking for the infrared 
asymptotic. We can separate the integral over fast variables and the density integral now 
becomes 



exp 



(2^ 



exp (In M) 



(35) 



where M is the integral over the fast variables, which in the first non-trivial order in per- 
turbation is a Gaussian integral without any prefactor 
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M 



V' 
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V 



exp 



(fk' 
(2^ 



cos2/3(0J(xi) + 0{(xi)) - cos2/?(0^(xi) +4{xi) 



COs2/3U'i(x2) +0{(X2)) -COS2/?(0^(X2) +0|(X2) 



(36) 



In functional integral Eq. (|35|) In M plays a role of renormalisation of the imaginary part of 
the kinetic operator. Now we calculate the functional integral and expand the result in the 
boson fields (p'f to get their bilinear combination, which gives the renormalisation of r]. We 
get 



M = ^^-^ J d'x,d^X2 



0a(a;i)7^a/3(a;i2)0^(x2) 



(37) 



where as usual xu = xi — X2- The elements of the matrix TZapix) are the exponentials 
of the boson Green's functions in the space-time representation. They are presented in 
Appendix B, where we give the detailed renormalisation calculation. The boson Green's 
functions entering the matrix TZaisi^u) are now different from those in Appendix A, due to 
the self-energy Eq. (|3lD. The retarded function in the momentum space is now 



Dr{uj, q) 



_ g2 _|_ iJfjiT,{uj, q) 



(38) 



and the corresponding functions DA{yj,q) and Dpioj.q) acquire the same self-energy. Note, 
that TZajsixu) is expressed via boson Green's functions in the not rotated Keldysh basis. 
However, the coefficient t] is easier to extract from the Green's functions in the rotated 



18 



Keldysh basis Eq. (|52D. The element TZr of the matrix 7laf3 in the rotated basis which 
corresponds to the renormahsation of the retarded function is 



= ^ ( 7^ll (x) - 7^22 (x) - (x) + 7^2l (x) ) . (39) 



If the cut-off is large, vpko > T, this yields after the Fourier transform 



TZr = iuj--^exp{-—\nkoa), (40) 



SO that the change in 77 is 



since /3 = V47r. Here as always, a is the lattice spacing which cuts-off the large momen- 
tum integration. So here we have a power-law rise in r^. When the cut-off is less that the 
temperature we get exponential renormahsation: 



T \ V 



Ari = riT exp ( 4—^ ) — — (42) 



Here tit is the effective damping at the scales vpko ~ T. That means, that on the scales of 
momenta less that the temperature, the imaginary part coefficient rj grows very rapidly, and 
we immediately get the diffusive dynamics. The renormahsation group treatment breaks 
down when the correction Eq. (^) becomes of the order of unity. That determines the 
length scales, where the dynamics becomes diffusive. The mean-free path is 
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(43) 



where rjo is the original value of the coefficient 7], proportional to the spin-phonon coupling 
constant (Eq. 

The estimate of the mean-free path allows us to estimate the diffusion coefficient as 
D = Ivp, so that 



We now briefly review our results. We started from ID Heisenberg Hamiltonian Eq. 
(0). Out goal was to calculate the spin-spin correlation function at long times and non- 
zero temperature, and check whether or not it had a diffusive pole Eq. (P in some region 
in phase space. To perform such a calculation we mapped our original problem onto ID 
fermion model Eq. (^) using the Jordan- Wigner transformation Eq. (^. We used the usual 
diagrammatic technique to calculate a few ffist orders of the perturbation theory (shown 
on Fig. 2 and Fig. 3). In order to go beyond that simple approximation we bosonised 
the fermion model. To get the long-time asymptotic of the spin-spin correlator, we had 
to combine bosonisation with Keldysh technique, which allows one to calculate real-time 
correlation functions at finite temperatures without any need in analytic continuation. 

The bosonisation procedure allowed us to find the general form of higher order corrections 
and to sum up the perturbation series. It turned out that in each order of perturbation the 




(44) 



VI. CONCLUSIONS 
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correction to the spin-spin correlator had the same form Eq. p6|), therefore we concluded, 
that the exact spin-spin correlation function do not acquire the diffusive pole from the 
summation of the perturbation series. This result is also known in the theory of the Sine- 
Gordon model (which coincides with our boson model), where it has been found that due 
to the infinite number of conserved charges excitations are propagating and not diffusive. 

We also check whether these results (the absence of spin diffusion in the perturbation 
theory) are robust with respect to small dissipation effects present in real physical systems. 
Specifically, we considered the effect of a weak spin-phonon interaction Eq. (p9D . We mapped 
the full problem (including the spin-phonon interaction) to the bosonic model. We found 
that the interaction with phonons leads to the boson self-energy Eq. (|3l|), the imaginary 
part of which is proportional to the constant, rj, at small momenta. Further, we applied the 
renormalisation group analysis; we integrated out "fast" degrees of freedom and showed that 
this constant f] grows moderately while the cut-off is larger than temperature (Eq. (|41|)), but 
grows exponentially after the cut-off becomes smaller than the temperature (Eq. (|4^)). At 
scales where the imaginary part of the spin-spin correlator becomes of the order of the real 
part the spin dynamics becomes purely diffusive. By associating the mean-free path with the 
scale on which the renormalisation procedure breaks down (namely, the renormalisation Eq. 
(^) becomes of the order of unity) we estimate the diffusion coefficient Eq. (0). Restoring 
the original units and estimating the spin-phonon coupling constant J' from the expansion 
of the exchange integral Eq. (pSj), we have 



where a is the lattice spacing and c is the speed of sound. 

Thus we found that the presence of the spin-phonon interaction changes the long-time 
behaviour of the spin-spin correlator, which becomes diffusive. 
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We present here the boson and fermion Green's functions, which we are using in our 
calculations. It is the certain similarity between them, that had inspired the bosonisation. 
As in the usual zero-temperature bosonisation we need the Green's functions in the space- 
time representations. We perform the bosonisation in the original Keldysh basis, but for 
simplicity we here calculate Green's functions in the rotated basis and them transform them 
back. 

The retarded fermion Green's function in the momentum space is 



where "-" is for the "left" and "+" for the "right" movers. For the simplification of the 
formulas the Fermi velocity is set equal to unity. The Fourier integral, which one has to 



large momenta. As usual in ID calculations we introduce a cut-off a by adding the exponent 



APPENDIX A 




(46) 



calculate in order to transform the function Eq. (^) to the real space is formally divergent at 



e "''^l to the integral. Thus we have 
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^ ' ^ ' (27r)' eTk + iO ^ ^ 



Now we have the perfectly converging integral and get 



GRix,t) = -^- , (48) 



27r (x T t)' + a2 

where 



eit) 

Similarly, the advanced function is 



1 , if t>Q 
, if t<0 



GaM = \^- ^ (49) 



The third Keldysh function in this basis for the right movers in momentum space is 

k 

F{e,k) = -2nita-nh—6{e- k). (50) 
The Fourier integral in converging and we get 



The function for the left movers has the opposite sign of the time variable. To be more 
careful with the pole one has to add ia to the space coordinate in Eq. 

That completes the calculation of the fermions Green's functions in the rotated Keldysh 
basis 

Gr 

Grot =1 I . (52) 

Ga F 
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To return to the original basis, which is needed for bosonisation, one has to perform the 
rotation 

/ G++ G+- \ 

G=\ \= RGrotR-\ (53) 

\G-+ G— I 

where the rotation matrix is given by 



1 



1 1' 



R = ^\ • (54) 



V2 



-1 1 



We now calculate the boson functions. We start again with the rotated basis. The 
retarded Green's function in the momentum space is 



DR{u;,q) = - l—r-. (55) 



Again we have to introduce the cut-off a. It is essential for the purposes of bosonisation to 
do it in exactly the same way as for the case of fermions, Eq. (^Tf). That way we get 



^ , , ^ X , ix — t + ia){x + 1 — ia) 

Dnix, t) = -— eit) In ) . . . 56 

4.1T [X — t — la) [x + t + ta) 



For the advanced function we get the same logarithm 



DA(x,t) = — ei-t) \n) -47 -4- (57) 

^ ' ' 47r ^ ^ {x-t-ia){x + t + ia) ^ ^ 
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The third function in this basis, Dp contains delta function as it's fermion counterpart Eq. 

m 

Df{u;, q) = -— coth ^ {6{u; - q) + 6{u; + q)) , (58) 



but we have to introduce the cut-off here. In real space we get 



Dp(x, t) = — In , ^, (59) 

' ^ 271 sinh7rT(x-t)sinh7rr(x + t)' ^ ^ 



where the zeros of the denominator should be treated in exactly the same way as in the 
fermion case ( see text after Eq. (^) ). 

One can clearly see, that the fermion functions and the arguments of the logarithm in 
the boson functions are constructed from the same elements. That is why the bosonisation 
works. To get exactly the same results in the fermion and boson perturbation series we 



have to turn to the original Keldysh basis Eq. (53) and calculate physical quantities, like the 



density-density correlation function, which is independent of Keldysh indices and therefore 
of the choice of the calculational technique. Then the boson and the fermion versions are 
the same in the limit a = (the physical quantities should not depend on that cut-off). 



APPENDIX B 

We start from the functional integral for the density-density correlation function Eq. 
(|ll|), where the kinetic term contains now non-zero imaginary part Eq. (selfen). We are 
interested in the infrared asymptotic of that correlator. Therefore we divide the boson field 
0^ in two part - "fast" and "slow" 

0^ = 0^ + (60) 

The integral Eq. (0) then becomes 
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nix, -X2) = - jV [0^ 0^] V 0{, (Pi 9i0^*(xi)r„^9i0^(x2) exp(-^) 



(61) 



where the prefactor contains only slow degrees of freedom since we are looking for the infrared 
asymptotic. We can separate the integral over fast variables as 



1 



V 



^14 



exp 



(2^ 



-r'(cos 2f3(f)l cos 2(3(j){ - sin 2(3(f)l sin 2(3(f)( - cos 2/?0^ cos 2/30^ 



+ sin 2/30^ sin 2/50^)]) . 



(62) 



The density integral becomes now 

nix,~X2) = ^ jv[<P{,ct>l] 9i0^*(xi)r„^ai0^(a:2) 



exp 



dPk' 
(2^ 



exp (In M) 



(63) 



where InM plays now a role of renormalisation of the imaginary part of the kinetic operator. 

We now calculate the integral over the fast degrees of freedom Eq. (|62D. We expand the 
exponent in Eq. (|6^) up to the first non-trivial order in perturbation (which is actually all we 



need, noting the results of Section |rV| ) and get the Gaussian integral without any prefactor 

<fk' 



y,2 

M = — I V 



exp 



{27ry 



d xid X2 



cos 2/50^ (xi) cos (xi) — sin 2/50^ (xi) sin (xi) — cos 2/902(xi) cos 2/502 (xi^ 
+ sin 2/50^(xi) sin 2/30^(xi)) 



cos2/30^(x2) cos2/30((x2) - sin 2/501 (X2) sin2/50((x2) - cos2/302(x2) cos 2/30^(^2) 



+ sin2/30^(x2)sin2/30^(x2)) 



(64) 
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We can now calculate the functional integral. The result is 



M ^ — J (Pxi(fx2< 



cos2/3(pl{xi) cos2/3(pl{x2) + sm2/3(pl{xi)sin2/3(pl{x2) 



exp if3^ / :.D, 

\ Jk'>ko (27r)' ' 



2 



1 - e 



ikx 



1-e 



-ikx^ 



+ 



COS 2/3^2(^1) COS 2/3^2(^2) + sin 2/302(^i)sm2/302(3^2) 
cos 2/901 (xi) cos 2/302(2^2) + sin 2/30^ (a;i)sm2/302(2^2) 



\ Jk'>ko [277) ^ 



-^ikx _ J)+-^-ikx 



COS 2/302 cos2/30i(a;2) + sin2/302(a:^i)-sm2/30i(a;2) 



where a; = — and is the momentum cut-off, delimiting fast variables from slow. 

The remaining integrals are similar to those, calculated in the regular perturbation se- 
ries. We have the Fourier transform of the boson Green's functions to the real space in the 
exponent, then take the exponential and Fourier transform back to the momentum space. 
The difference is that we now have another set of Green's functions - with non-zero imagi- 
nary part - and the momentum integration in limited by the cut-off fco- The result of this 
integration will now depend on the cut-off. For the retarded and advanced functions we get 



Dn{x,t-M) = ^ m e-'^''2iImlE, 



ko {rjt — i{x + t)) 



ko {rjt + i{x — t)) 



(66) 
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I 

An 



-t) e-'^'^'hilml El 



ko {r]t - i{x + t)) 



ko {r]t + i{x - t)) 



(67) 



where Ei is the exponential integraL The result for the third Keldysh function Dp depends 
on whether the cut-off is larger or smaller than the temperature. For ko > T we get 



i 1 

Df{x, t;ko> T) = -— -Re 

ZTT 2 



El 



ko {a + i{x + t)) 



+ Ei 



ko (a + i{x — t)) 



(68) 



so that at the origin 



Df{0; ko>T) = — ln(fco«) 



(69) 



which is a small number. For the other case, ko < T, in the limit of large x and t and 
assuming that original rj is much less than the temperature we have 

iT ( sin ko{x + t) sin ko{x-t)\ 

and at the origin is 

DFiO;ko<T) = --^, (71) 
nko 

which is extremely big. 

We can now proceed with the renormalisation of the imaginary part coefficient r/. To do 
that we expand Eq. ( p5D in boson fields and get their bilinear combination 
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0i(xi)0i(x2) exp [ip^ 



k'>ko (2n 



2 



DZ- ( 1 - e^^^n ( 1 - e-'^"" 



+02(3^1)02(^2) exp ij3 



k'>ko {2tt 



,2 k 



2 1 - - D++ d-e''"') 



\ Jk'>ko [ZTCj ^ 



6^(xi)0^(x2)exp {i(3 



k'>kQ {2tt) 



(72) 



which can be abbreviated as 



M = J d'xid'x2 



0^(a;l)7^^/3(a;)0^(x2) 



(73) 



where as usual x = xi — X2- 

The imaginary part coefficient rj is most transparent in the rotated Keldysh basis Eq. 
(^). The element of the matrix TZajsix) in the rotated basis which corresponds to the 
renormalisation of rj in the retarded function in terms of original elements is 



= 2 ( ^ii(^) - ^22(3;) - 7^l2(x) + 7^2l(a;) 



(74) 



For the case of the large cut-off this yields after the Fourier transform iuji]-^ exp(— ^ In k^a), 
so that the change in 77 is 
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(75) 



since (3 — \f^. Here as always, a is the lattice spacing which cuts-off the large momen- 
tum integration. So here we have a power-law rise in 77. When the cut-off is less that the 
temperature we gain different exponential, so that now 

A, = , exp (4I] (76) 

That means, that on the scales of momenta less that the temperature, the imaginary part 
coefficient?] experiences a tremendous growth, and we immediately get the diffusive dynam- 
ics. 
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